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Effects of bedrock lithology and subglacial till on the motion of Ruth Glacier, Alaska, 

deduced from five pulses from 1973-2012 

ABSTRACT 

A pulse is a type of unstable glacier flow intermediate between normal flow and surging. Using 
Landsat MSS, TM, and ETM+ imagery and feature tracking software, a time-series of mostly 
annual velocity maps from 1973 to 2012 was produced that reveals five pulses of Ruth Glacier, 
Alaska. Peaks in ice velocity were found in 1981, 1989, 1997, 2003, and 2010; approximately 
every 7 years. During these peak years the ice velocity increased 300%, from approximately 40 
m/yr to 160 m/yr, and occurred in an area of the glacier underlain by sedimentary bedrock. Based 
on the spatio-temporal behavior of Ruth Glacier during the pulse cycles, we suggest the pulses 
are due to enhanced basal motion via deformation of a subglacial till. The cyclical nature of the 
pulses is theorized to be due to a thin till, with low permeability, that causes incomplete drainage 
of the till between the pulses, followed by eventual recharge and dilation of the till. These 
findings suggest care is needed when attempting to correlate changes in regional climate with 
decadal-scale changes in velocity, because in some instances basal conditions may have a greater 
influence on ice dynamics than climate. 


INTRODUCTION 
Ruth Glacier 

Ruth Glacier is a 58 km long alpine glacier in the Alaska Range of central Alaska, USA 
(Figure 1), residing in a subarctic continental climate (Shulski and Mogil, 2009) with a 
predominantly southern aspect. Ruth Glacier has a multi-lobed accumulation area consisting of a 
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West Fork, Northwest Fork, and North Fork. The headwall of the Northwest Fork reaches 4300 
meters above sea level and is approximately 4 km from the summit of Mount McKinley, the 
highest peak in North America. The separate lobes of the accumulation area converge at the head 
of the Great Gorge, an approximately 12 km long, 1.5-2 km wide valley with 1500 m high steep- 
sided walls that forms the upper ablation area between the 12 km and 26 km marks (Fig. 1). In 
1983 the ice depth in the Great Gorge was measured using seismic methods and found to be up 
to 1150 meters thick (unpublished data of K. Echelmeyer). The lower ablation area, below the 
Great Gorge, is approximately 3 1 km long and 3-4 km wide and ends in a small piedmont lobe 
roughly 8 km wide, which is partially covered by vegetation. The average ice surface slope 
within the Great Gorge is 2.6°, and below the gorge the average ice surface slope is half as much, 
1.3°, per the National Elevation Dataset (Gesch and others, 2009). At the base of the Great Gorge 
there is some crevassing where the ice exits the gorge; the rest of the ablation area is free of 
major crevassing (as seen in the 30-meter Landsat imagery used in this study), but the medial 
moraines exhibit a slightly wavy pattern. The Great Gorge, and parts of the West Fork, 
Northwest Fork, and lower ablation area are underlain by Tertiary (Paleocene) biotite -muscovite 
granite and quartz monzonite of the McKinley Sequence (delineated by the red polygons and 
denoted by the symbol Tpgr in Figure 1) (Reed and Nelson, 1980; Gamble and Reed, 1996). The 
North Fork and the ablation area immediately below the Great Gorge are underlain by 
sedimentary rocks from the Cretaceous and Jurassic periods (denoted as KJf in Figure 1). These 
rocks are part of the Kahiltna Flysch Sequence and are comprised of argillite, fine to coarse 
greywhacke, conglomerates, and thin layers of chert and limestone (Reed and Nelson, 1980; 
Csejtey and others, 1992). Part of the piedmont lobe is underlain by Tertiary (Miocene) 
sedimentary bedrock comprised of sandstone, siltstone, shale, and claystone of the Tyonek 
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Formation (delineated by the yellow dotted polygon and denoted as Tty in Figure 1), which is 
part of the Kenai Group (Solie and others, 1991). 

Little has been published concerning Ruth Glacier. Mayo (1978) notes Ruth Glacier is a 
pulsing glacier and mentions the existence of small potholes in its surface that are the remnants 
of large crevasses. Hall and Ormsby (1983) examined Ruth Glacier using Landsat Multispectral 
Scanner and SEAS AT synthetic aperture radar data from the summer of 1978 and noted the fim 
zone was located within the Great Gorge. Heid and Kaab (2012a) produced velocity 
measurements for nine glaciers in the Alaska Range using optical feature tracking methods and 
noted only Ruth Glacier had accelerated during the period from 1986/87 to 2009/10; the zone of 
acceleration was restricted to the ablation zone below the Great Gorge. Burgess and others 
(2013) used synthetic aperture radar offset tracking to measure the wintertime speed of all 
glaciers in the Alaska Range during the late 2000s; the exact year is not given for each glacier, 
varying from 2007 to 2010. It was found that Ruth Glacier was moving at 0.6 to 1.0 m/day (219 
to 365 m/yr) in the Great Gorge, with the ice gradually decelerating from approximately 0.5 
m/day (182.5 m/yr) immediately below the Great Gorge to being stagnant at the terminus. Ward 
and others (2012) note the igneous bedrock (Tpgr) of the Denali massif is significantly less 
fractured than the surrounding tectonically deformed sedimentary bedrock (KJf and Tty). 
Fracture spacing in the granitic bedrock is on the order of decimeters to hectometers, while 
fracture spacing in the sedimentary bedrock is on the order of centimeters to decimeters. This 
difference in fracture spacing causes a large difference in rates of glacier erosion (plucking) 
between the two predominant bedrock types. The result is that the valley in which Ruth Glacier 
resides tends to be narrow, with steep-sided walls (often 60° or greater in slope) in regions 
underlain by the granitic bedrock (eg. The Great Gorge and between the 41 and 48 km marks), 
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and the valley tends to be wider with shallow-sloped walls (often 45° or less) when underlain by 
sedimentary bedrock (Ward and others, 2012). Thus, the geometry of the glacier, and 
consequently ice thickness and basal shear stress, are dictated partly by the underlying bedrock 
type. 

Pulsing glaciers 

Mayo (1978) defines glacier pulses as periodic unstable flow that is lesser in magnitude 
than surges, and therefore, pulsing glaciers are intermediate between nonnal and surging 
glaciers. Pulsing glaciers may be identified by the existence of characteristic wavy medial 
moraines, large-scale wavy foliation, or boudinage that may be present on only part of the 
glacier, indicating only part of it pulses. If the characteristic pattern is continuous, it is suggested 
the glacier experiences regular, periodic pulses. Mayo (1978) notes that potholes (remnants of 
crevasses) probably form in the zone of ice loss on pulsing glaciers. Approximately 140 glaciers 
in Alaska have been identified with these characteristics in aerial photography from the 1960s 
and 1970s by Austin Post and L.R. Mayo. Most of these pulse-type glaciers occur in the same 
regions as surge-type glaciers, in the Alaska Range, Chugach Mtns., St. Elias Mtns., and 
Wrangell Mtns. Pulsing behavior may occur in conjunction with the drainage of large glacier- 
dammed lakes, suggesting a link between basal hydrology and pulsing flow. Mayo (1978) lists 
some pulsing glaciers and the dates of their observed pulses: West Gakona Glacier (1949), 
Nizina Glacier (1961), Trident Glacier (1970, 1971), MacLaren Glacier (1971, 1972), and 
Tokositna Glacier (1971, 1972). Additionally, undated pulses are listed for: Capps Glacier, 
Eureka Glacier, Hayes Glacier, Kahiltna Glacier, Miles Glacier, Ruth Glacier, Sanford Glacier, 
and the West Branch of Sheridan Glacier. 
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Sometime between 1974 and 1977 Trapridge Glacier, Yukon Territory, Canada, began a 
prolonged surge that lasted until 2005 (Frappe and Clarke, 2007). During the surge a string of 
five 4-year pulses were detected from 1981 to 2002, with amplitudes of roughly 10 m/yr, a 33% 
increase in velocity. Trapridge Glacier is a polythermal glacier with temperate upper ice and 
subzero basal ice, and it is underlain by a deformable, permeable till up to 10 m thick (Blake and 
others, 1992; Stone, 1993). Increased ice velocity during the surge, and by inference during the 
pulses, is attributed to sliding at the bed caused by failure of the basal till. Till failure reduces 
basal friction and transfers stress laterally to the glacier margins, thereby producing an area of 
plug-flow along the central portion of the glacier, allowing ice there to accelerate. During a 
similar time frame, from approximately 1980 to 2002, Black Rapids Glacier, a temperate glacier 
in the central Alaska Range (Harrison and others, 1975), experienced two consecutive 12-year 
pulses, with amplitudes of roughly 20-25 m/yr, representing velocity increases of 55-65% 
(Nolan, 2003). Black Rapids Glacier is underlain by a till 5-7 m thick (Nolan and Echelmeyer, 
1999), and it has been shown via borehole tiltmeter measurements that up to 70% of the ice 
velocity during the pulses was due to deformation within the till at depths >2 m below the ice/till 
interface (Truffer and others, 2000). Again, it was theorized that till failure reduced stresses in 
the center of the glacier and transferred them towards the margins. This process promotes 
increased basal motion along the glacier centerline and increased ice deformation near the 
margins (Truffer and others, 2001; Nolan, 2003). 

Ice dynamics influence the length, area, and volume of a glacier. Changes in ice 
dynamics, therefore, necessarily result in changes in the geometry of a glacier, which in turn 
affects its mass balance. For example, pulses and surges may rapidly transfer ice to lower 
altitudes where it is subject to increased air temperature and ablation, thereby promoting a 


5 


115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 


reduction in mass balance. Knowing that dynamic instabilities such as pulses and surges are 
intimately tied to the basal conditions of a glacier, an improved understanding of the relationship 
between basal conditions and ice dynamics will further our understanding of the influences on a 
glacier’s mass balance. In this study we will examine the cyclical dynamic behavior of Ruth 
Glacier to determine the amplitude, frequency, and spatial extent of its pulses. Based on the 
spatial and temporal structure of the velocity fields during multiple pulse cycles, we will infer the 
basal conditions beneath Ruth Glacier and describe the influence of bedrock lithology on the 
behavior of its pulses. 


METHODS 

Optical feature tracking 

Optical feature tracking methods were used to produce a time-series of ice surface 
velocity maps spanning 1973-2012 for Ruth Glacier in order to delineate its periodic dynamic 
behavior. COSI-Corr feature tracking software (Leprince and others, 2007) was used in 
conjunction with Landsat Multispectral Scanner System (MSS), Thematic Mapper (TM), and 
Enhanced Thematic Mapper Plus (ETM+) imagery to produce the velocity maps (see Table 1 for 
image dates). COSI-Corr was chosen due to its proven precision (Scherler and others, 2008; Heid 
and Kaab, 2012b), its ability to produce accurate displacement measurements in areas of low 
contrast or with light cloud cover (Heid and Kaab, 2012b), and its ease of use. COSI-Corr 
produces an intitial estimate of the displacement of features between sequential images by use of 
the Fourier Shift Theorem which states the displacement is found in the phase portion of the 
normalized cross power spectrum computed from the Fourier transform of each image 
(Shekarforoush and others, 1996). Equal-sized subsets from each image centered on the same 
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pixel, called reference and search windows, corresponding to the temporally first and second 
images, respectively, are matched to one another and a displacement measurement produced with 
1 -pixel accuracy. A re- weighted least squares phase minimization algorithm is then applied to 
determine the final displacement. This algorithm minimizes the phase difference between the 
initial displacement estimate and an ideal theoretical one, thus producing subpixel accuracy. This 
process is repeated in a grid-like pattern across the entire image, producing individual 
displacement measurements at predefined intervals. The displacement measurements (vectors) 
produced by COSI-Corr are post-processed using a neighborhood analysis routine to remove 
erroneous vectors. Any vector that is more than ±2 standard deviations away from the mean 
magnitude or direction of its local neighborhood, consisting of at least nine adjacent vectors, is 
deemed an anomaly and removed. Afterwards, the vector field is visually inspected and any 
remaining anomalies are manually removed. The individual vectors are interpolated using a 
linear inverse-distance weighting scheme to produce a velocity raster. In May, 2003, the scan- 
line correction mirror of Landsat-7 ETM+ failed, resulting in data voids in images acquired 
afterwards; these images are commonly called SLC-off images. Fortunately, the portion of Ruth 
Glacier in which the pulses are tracked (the ablation zone below the Great Gorge) lies near the 
center of the Landsat-7 images where the data voids do not exist. Therefore, more advanced 
feature tracking techniques, such as Ahn and Howat’s (2011) null exclusion method or 
orientation correlation (Fitch and others, 2002; Haug and others, 2010) were not necessary when 
using SLC-off images for Ruth Glacier. 

Accuracy 

The accuracy of the velocity fields produced from feature tracking algorithms depends 
on the precision of the matching method and how well the two images being matched are aligned 
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to one another. Heid and Kaab (2012b) evaluated the precision of COSI-Corr using Landsat-7 
ETM+ 15-meter panchromatic imagery and found the root mean square error (RMSE) of its 
displacement measurements to be ±0.9 meters in the x-direction and ±0.8 meters in the y- 
direction, for a total RMSE of ±1.2 meters, or ±0.08 pixels. If we apply this same error (±0.08 
pixels) to 30-meter imagery, the expected error is ±2.4 meters, and for 60-meter imagery the 
expected error is ±4.8 meters. Storey and Choate (2004) have shown the RMSE of the geometric 
accuracy of Landsat-5 TM data is ±5.5 meters (±0.2 pixels), and Lee and others (2004) showed 
the average geometric error of Landsat-7 ETM± imagery is less than ±5 meters (±0.166 pixels). 
Table 2 (last column) shows the listed geometric accuracy for each Landsat image used in this 
study, upon receipt from the U.S. Geological Survey; the mean of the 22 values for TM/ETM± 
images is 0.16 ± 0.04 pixels (4.8 ±1.2 meters) (mean ± 1 standard deviation), which compares 
well with the values given by Storey and Choate (2004) and Lee and others (2004). For the 12 
MSS images, the mean geometric error is 0.47 ± 0.12 pixels (28.2 ± 7.2 meters). 

For each image pair vectors on dry, snow-free, cloud-free land were analyzed to 
determine the mean georeferencing error between the two images (Table 2). The average of these 
mean georeferencing errors is 0.22 ± 0.09 pixels (7.0 ±3.0 meters) for 30-meter TM and ETM± 
imagery and 0.45 ±0.15 pixels (24.1 ± 6.7 meters) for MSS imagery. The mean georeferencing 
error for each image pair was subtracted from the on-ice vectors to improve the accuracy of the 
velocity fields. The mean georeferencing error was also subtracted from the off-ice vectors, 
resulting in a residual error between images which represents their final misalignment (Table 2). 
We note the average of the mean residual errors for the TM/ETM± image pairs, 0.15 ± 0.07 
pixels (4.9 ± 2.6 meters), compares well with values quoted above by Storey and Choate (2004), 
Lee and others (2004), and with the mean geometric accuracy for an individual image. Thus, 
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after removal of the mean georeferencing error, the TM/ETM+ image pairs are aligned to one 
another as precisely as each individual image is aligned to the ground control points used to 
geolocate it. We also note the image-to-image misalignment for each TM/ETM+ image pair is 
within one standard deviation of the precision of COSI-Corr (±0.08 pixels), which suggests that 
after removal of the mean georeferencing error the residual misalignment of the TM and ETM+ 
image pairs approaches the limits of detection by COSI-Corr. The mean georeferencing error for 
MSS data, as measured by COSI-Corr (Table 2) is 0.45 ± 0.15 pixels (24.1 ± 6.7 meters), and 
after removal of this error, the residual image-to-image misalignment is 0.25 ± 0.07 pixels (14.1 
± 4.2 meters). Thus, after removal of the georeferencing error, the MSS image pairs are aligned 
more precisely to one another than each individual image was originally aligned to its ground 
control points upon receipt from the U.S. Geological Survey (compare the average residual error 
of 0.25 pixels to the average geolocation error of 0.47 pixels). Using the root sum of squares 
method, we estimate the combined error in the velocity fields, resulting from the mean residual 
error and COSI-Corr’s precision, to be 5.1 meters (0.17 pixels) for the 30-meter Landsat 
TM/ETM+ data, and 15.7 meters (0.26 pixels) for the 60-meter Landsat MSS data. 


RESULTS 

COSI-Corr was unsuccessful in producing reliable displacement measurements within the 
Great Gorge and in the accumulation area due to snow cover and a lack of surface features to 
track, but it was quite successful below the Great Gorge. A time-series of 29 velocity maps 
(Figures 2 to 5) shows the evolution of the pulses of Ruth Glacier in time and space, below the 
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Great Gorge. From 1973-74 to 1976-77 (Fig. 2a-c) much of the ablation area below the Great 
Gorge is moving at 75 m/yr or less. Then in 1977-78 (Fig. 2d) the ice down to the 30 km mark 
accelerates to almost 200 m/yr. In 1978-80 (Fig. 2e) the ice between the 35 km and 40 km marks 
has accelerated relative to previous years, and by 1980-81 (Fig. 21) all the ice down to the 45 km 
mark has accelerated, much of it moving 150 m/yr or more, indicating a doubling of ice velocity 
compared to the period from 1973 to 1977. From 1981-1982 (Fig. 2g) to 1983-84 (Fig. 3a), the 
glacier gradually decelerates and returns to its pre -pulse velocity. In Figure 3 there are two pulses 
evident between 1984-85 and 1995-99 that show a similar pattern of gradual acceleration and 
deceleration of the ice, with peak velocities in 1987-91 (Fig. 3e) and 1995-99 (Fig. 3h). Due to 
lack of cloud- free imagery between 1987 and 1991, and lack of data acquisition during the late 
1990s, these pulses are not well defined in time, but it is clear the two pulses did occur due to the 
increased velocities during these four-year time spans. From 1999-00 to 2001-02 (Fig. 4) the 
cyclical pattern has continued, with a gradual increase in velocity from roughly 50-75 m/yr in 
1999-00 (Fig. 4a) to >100 m/yr in 2001-02 (Fig. 4c) between the 26 and 35 km marks. By 2002- 
03 (Fig. 4d) velocity reached a peak, with characteristic values of 150-200 m/yr between the 26 
and 42 km marks, indicating the ice has more than doubled its speed compared to three years 
earlier. Afterwards, the ice decelerates to its pre-pulse velocity by 2005-06 (Fig. 4g). A fifth 
pulse is evident in Figure 5, which shows an acceleration of the ice from 2006-07 (Fig. 5a) to 
2009-10 (Fig. 5d), with the peak being in 2009-10, followed by deceleration in 2010-11 and 
2011-12 (Fig. 5e,f). Characteristic peak ice velocities during this last pulse, between the 26 and 
42 km marks, are 150-200 m/yr, indicating a three- to four-fold increase in velocity compared to 
2006-07. 
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To help illustrate the evolution of the pulses through time, the velocity at two points 
along the glacier, 36 km and 43 km, are plotted in Figure 6. This graph clearly illustrates the 
periodic nature of the pulses and highlights the five individual peaks. It should be noted that for 
plotting purposes, velocity was assigned to the 2 nd year of each image pair used to produce a 
velocity field, or in the case of the 1987-91 and 1995-99 velocity fields, velocity was assigned to 
a middle year, 1989 and 1997, respectively. Inspection of this graph yields characteristic pulse 
amplitudes (peak value minus trough value) in excess of 120 m/yr (a 300% increase in velocity) 
at the 36 km mark, and 60-70 m/yr (a 50-75% increase in velocity) at the 43 km mark for the 
1981, 2003, and 2010 peaks. The amplitudes of the 1989 and 1997 pulses are not as great 
because their respective velocity fields were produced using image pairs spanning four years, 
thus the peak velocity values are averaged with surrounding years resulting in lower amplitudes. 
Based on the occurrence of the peaks as described here, the pulses of Ruth Glacier have a 
characteristic period of approximately 7 years. 

To further aid in analyzing the pulses, longitudinal profiles of velocity, acceleration, and 
strain rate were produced for each velocity field (Figures 7-10). The longitudinal velocity and 
acceleration profiles in Figures 7a and 7b clearly show the extent of the pulse from 2009 to 201 1 . 
All the ice below the Great Gorge down to roughly the 52 km mark noticeably accelerated during 
this time. It is worthwhile to note that in the area of the 27 km mark, the 40-42 km marks, and 
the 49 km mark, there are inflections in velocity and acceleration, and prominent peaks (or 
troughs) in strain rate in 2009-10, 2010-1 1, and 2011-12. The location of these inflections, peaks, 
and troughs approximately coincide with the transition from granitic to sedimentary bedrock at 
the 28 km mark, a change from sedimentary to granitic bedrock at the 41 km mark, and a change 
back to sedimentary bedrock at the 48 km mark, as seen in Fig 1. Examination of the 


11 


252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 


acceleration profiles (Fig. 7b) reveals the ice from the 26 km mark to roughly the 50 km mark 
accelerates (in varying amounts) simultaneously; there is no evidence of a wave front 
propagating downglacier. The strain rate between the 37 and 43 km marks changes from 
predominantly compressive before the pulse, to tensile during the pulse, and back to compressive 
afterwards; again, indicating the ice in this area is responding simultaneously to applied stresses. 
The longitudinal profiles in Figure 8 span the years 1999 to 2006 and highlight the behavior of 
Ruth Glacier immediately before, during, and after the 2002-2003 pulse. As with the 2009-2010 
pulse, there are noticeable inflections in velocity and acceleration, and noticeable peaks or 
troughs in strain rate, near the transitions between granitic and sedimentary bedrock at 28 km, 41 
km, and 48 km. Again, it is evident that the entire stretch of ice between 26 and 50 km has 
accelerated simultaneously, and the same pattern of changing strain rate, from compressive to 
tensile, and back to compressive, is evident before, during, and after the pulse, respectively. The 
2003-04 acceleration profile (Fig. 8b) shows a well-defined area of acceleration between the 42 
and 50 km marks, atop granitic bedrock; this same area also increased in acceleration in 2010-1 1 
(Fig. 7b) while the rest of the ice upglacier decelerated. Both of these acceleration events 
occurred one year after, and downglacier, of the peak pulse events in 2003 and 2010, which 
occurred atop sedimentary bedrock. Figure 9 shows the longitudinal profiles for the years 
surrounding, and including, the 1987-91 and 1995-99 pulses. Although these two pulses are not 
as well defined as the 2002-03 and 2009-10 pulses, the same patterns are still evident. There are 
inflections in the velocity and acceleration profiles (Fig. 9a, b), and peaks or troughs in the strain 
rate (Fig. 9c), near the same locations as the later pulses, at 26 km, 40-42 km, and 50 km. Again, 
large portions of the ice between the 26 and 50 km marks accelerate and decelerate in unison. 
During the 1987-91 and 1995-99 pulses, there is no obvious acceleration of the ice between the 


12 


275 

276 

277 

278 

279 

280 

281 

282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 


42 and 50 km marks (Fig. 9b) (as seen in the later pulses (Fig. 7b, 8b)), probably due to the poor 
temporal resolution of these two pulses. Lastly, Figure 10 shows the longitudinal profiles for the 
years surrounding, and including, the 1980-81 pulse. Again, inflection points in the velocity and 
acceleration profiles are seen at, or near, the 26 km, 40-42 km, and 50 km marks (Fig. 10 a,b), 
but the pattern of peaks and troughs in strain rate seen in the other pulses is not clear. 


DISCUSSION 

To place our results in context with other surging and pulsing glaciers, we briefly discuss 
till deformation and flow instabilities with regards to glacier dynamics. In turn, these discussions 
help identify the likely basal conditions of Ruth Glacier during its pulses. 

Glacier motion by till deformation 

The location of surging glaciers in Svalbard is positively correlated with areas of fine- 
grained sedimentary bedrock, which are easily eroded (Jiskoot and others, 2000), while in NW 
North America surging glaciers are often found in fault- shattered valleys (Post, 1969). Glaciers 
atop these areas are more likely to develop a soft bed than when atop intact crystalline 
metamorphic or igneous bedrock. Sedimentary bedrock will produce a fine-grained till that will 
have a low penneability which may be more prone to instabilities (Murray and Porter, 2001). 
Instability within a till may develop due to decreased effective pressure, via increased basal 
water input. As long as the discharge of water at the glacier bed is small, steady-state water 
pressure will approximate ice overburden pressure and the water will flow by Darcian transport 
through the till, or by a unifonn water film at the base of the ice, and the till will remain 
consolidated. When pore -water pressure in the till equals or exceeds ice overburden pressure the 
till weakens and may dilate (Willis, 1995). Dilation of the till causes porosity and penneability to 
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increase, thus decreasing water pressure and shear strength (Murray, 1997). When a till dilates it 
can no longer support shear stress and it deforms, allowing the glacier to accelerate. The 
increased ice motion may occur as slip at the ice/till interface in which shear strain in the till 
occurs in the uppermost layer of sediment. Shear strain may also occur within the till along a 
plane of weakness, or it may be distributed throughout the till by the process of dilatant 
hardening. When a till dilates, pore volume increases, allowing water to flow towards the 
opening pores in the dilating sediment. If the rate of dilation is greater than the rate of pore-water 
flow then water pressure will decrease in the dilated sediment and it will harden (Reynolds, 
1885). When the shear band hardens, the strain is focused elsewhere within the till; thus, the 
strain is distributed from one plane of weakness to the next. 

Truffer and others (2000) suggest surges could be triggered by widespread dilation of 
basal till as shear stresses exceed a critical value, along with attainment of a glacier geometry 
that produces large basal shear stresses. Failed tills are concentrated under thicker ice, such as 
along the glacier centerline, where driving stress is greatest and may exceed till yield strength. 
Nolan (2003), when investigating the behavior of Black Rapid Glacier, proposed the idea of 
waves of till failure (activation waves) that propagate underneath a glacier as follows: 1) Failed 
till along a section of glacier centerline causes increased stress laterally and longitudinally 
because the driving stress is redistributed elsewhere; 2) The glacier speed increases locally due to 
failed till; 3) When the driving stress is distributed elsewhere, any till that was at, or near, a point 
of failure will fail, thus increasing the area of the glacier with increased speed; 4) This process 
creates a wave of till failure and subsequent ice motion that propagates downglacier via a 
positive feedback loop by continuously placing increased stress on the downglacier till; 5) As the 
area of failed till increases longitudinally, the glacier moves faster everywhere because of 
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reduced bed friction and greater basal motion. A similar mechanism was discussed by Frappe 
and Clarke (2007) to explain the surge of Trapridge Glacier. The till defonnation paradigm 
allows glacier velocity to increase greatly without requiring large changes in glacier geometry or 
the subglacial hydraulic system, assuming the till is near its failure threshold (Truffer and others, 
2001 ). 

Fowler and others (2001) noted that if the till layer has a low transmissivity (i.e. the till is 
thin) and low permeability then oscillations in effective pressure, ice thickness, and ice flux can 
occur. Ice speed during a pulse depends primarily on basal roughness. If the basal roughness is 
large, then ice flow is sufficiently restricted to allow oscillatory behavior to occur; otherwise, a 
major surge can occur if the basal roughness is small, because there is little restriction to ice 
flow. Fowler and others (2001) also suggest that if a glacier is temperate and resides atop a thin 
till, the ice may activate (i.e. the till fails beneath it) before any significant ice motion because 
the activation waves are faster than the ice motion; so, the ice slumps forward rapidly moving as 
a block, or plug. Therefore, there is no surge front that would otherwise separate fast-moving 
surging ice and slow-moving quiescent ice. Plug flow has been observed on Black Rapids 
Glacier (Harrison and Post, 2003) and Trapridge Glacier (Frappe and Clarke, 2007), both of 
which have surged atop a till. Once a surge has been initiated by till failure, less stress is required 
to keep the till in a failed state, allowing the surge to continue. This occurs because drag at the 
glacier base increases longitudinal shear stress that might cause a till to fail even when water 
pressures are less than the critical Coulomb threshold (Boulton and others, 2001). The 
accelerated ice motion will cease when the ice can no longer be shoved forward by driving 
stresses, either by obstacles restricting down glacier movement or by decreased upglacier 
stresses, then the basal shear stress decreases below the till failure point and the till heals. Or, the 
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accelerated ice motion will cease if drainage of water from the till occurs, reducing water 
pressure and allowing the till to heal. However, a till with low permeability may not drain 
completely; thus, the till may not fully heal and it is left primed for another pulse episode when 
enough water is introduced again (Nolan, 2003). A thin till, which has lesser water storage 
capacity than a thick till, will refill quickly and surpass its Coulomb failure threshold sooner than 
a thick till would, resulting in more frequent pulses. 

Basal conditions beneath Ruth Glacier 

Knowing that the sedimentary bedrock beneath Ruth Glacier is highly fractured and 
easily eroded, and the fact that Ruth Glacier exhibits a dynamic behavior similar to other glaciers 
kn own to reside atop deformable tills by accelerating and decelerating in a cyclical, plug-like 
fashion, it seems likely a deformable till exists beneath Ruth Glacier between the 28 km and 41 
km marks. The area from the 41 km mark to the 48 km mark is underlain by granitic bedrock that 
will not as readily erode to produce a till, but there is likely some sediment beneath the glacier in 
this area that was transported from above, but which is spatially sparse and does not promote 
widespread increased ice motion via till dilation and deformation. The fact that the ice in this 
area (lying atop the granitic bedrock) has accelerated during the pulses, although to a lesser 
degree than ice atop sedimentary bedrock, is likely due to stress transfer and longitudinal 
coupling with the ice upglacier (note the lesser velocity at the 43 km mark in Fig. 6 compared to 
the 36 km mark). The inferred existence of a deformable till requires the basal ice and the till to 
be at the pressure-melting point; therefore, we conclude Ruth Glacier is warmed-based below the 
Great Gorge. The simultaneous acceleration of the ice from the 26 km mark, immediately below 
the Great Gorge, to roughly the 50 km mark, and the lack of a wave front, provide further 
evidence that Ruth Glacier is wann-based in this area, because Fowler and others (2001) theorize 
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temperate glaciers atop deforming tills will slump forward, thus no wave front forms. The lack of 
a wave front is in contrast to observations of typical Alaska-style surging glaciers wherein a 
surge front separates quiescent ice from surging ice, as noted on Variegated Glacier (Kamb and 
others, 1985) and Bering Glacier (Roush and others, 2003; Turrin and others, 2013). The 
dynamic oscillating (pulsing) behavior of Ruth Glacier suggests the till is thin with a low 
permeability (Fowler and others, 2001). 

Distinguishing pulses from surges 

Alaska-style surges are controlled by the morphological evolution of the basal drainage 
system beneath temperate ice from a channelized system to a distributed cavity system (Kamb 
and others, 1985) and are traditionally assumed to occur atop a hard bed. Svalbard-style surges 
are controlled by the thermal evolution of the basal ice from sub-freezing to temperate (Fowler 
and others, 2001; Murray and others, 2003) and occur atop a soft bed (a basal till). Both Alaska- 
style and Svalbard-style surges can result in a surge front (kinematic wave) moving downglacier 
that separates quiescent ice, moving primarily by internal deformation of the ice column, from 
surging ice that moves primarily by slip at the ice/bed interface, or by deformation within a basal 
till (in the case of a Svalbard-style surge). Pulses are controlled by the evolution of the basal till 
from a healed state to a dilated state, and occur where the base is at the pressure-melting point. If 
the till beneath Ruth Glacier froze between pulses, then a surge front would form in the manner 
of a typical Svalbard-style surge during each pulse, but the lack of any wave front suggests the 
base is perennially warm. As noted above, Trapridge Glacier experienced five low-amplitude 
pulses during a prolonged surge (Frappe and Clarke, 2007), and Black Rapids Glacier, also a 
surge-type glacier, experienced two low-amplitude pulses during quiescence (Nolan, 2003). 
These studies suggest pulsing is a dynamic behavior apart from typical surging that may occur 


17 


390 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 


either during an active, low-amplitude surge or during quiescence. Ruth Glacier is not known to 
surge in the traditional sense, so pulsing may also occur on glaciers other than surge-type that 
have the requisite warm base and till, and adequate shear stress and basal water pressure to dilate 
the till. 

Implications for decadal-scale velocity measurements 

As mentioned above, Heid and Kaab (2012a) noted Ruth Glacier accelerated between 
1986/87 and 2009/10. This perceived acceleration is a result of the years chosen for use in their 
velocity determinations. The 1986/87 period falls between pulses, while the 2009/10 period is at 
the peak of a pulse (Fig. 3), thus there is a perceived acceleration of >120 m/yr during this 
period. If one had instead chosen images acquired in 1980/81 and 2007/2008 and perfonned the 
same analysis, the result would be a deceleration of >-120 m/yr, the opposite of what Heid and 
Kaab (2012a) reported. In either case, the pulses which occurred during the interval between 
velocity measurements are missed. Optical feature tracking methods have now matured to the 
point where decadal-scale regional, and perhaps global, studies of ice dynamics are possible, as 
exemplified by Heid and Kaab (2012a). When performing such work, it is important to be aware 
of any surging or pulsing glaciers in the chosen study area and how their periodic behavior might 
affect results. In some instances, such as on Ruth Glacier, basal conditions may exert a greater 
influence on ice dynamics than regional changes in climate on annual and decadal time frames. 

CONCLUSIONS 

From 1973 to 2012, a span of 39 years, Ruth Glacier experienced five pulses that have 
gone unnoticed until now. These pulses are of low amplitude compared to typical surges of 
glaciers in Alaska, and they occur without causing an advance of the terminus and without 
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widespread crevassing, perhaps helping to explain their oversight. The pulses have occurred on a 
regular basis, approximately every seven years, with peaks in 1981, 1989, 1997, 2003, and 2010, 
and with increases in ice velocity of >120 m/yr during the peaks in 1981, 2003, and 2010. The 
pulses occur in an area of the glacier underlain by sedimentary bedrock, and we suggest that the 
pulses are a result of enhanced basal motion due to deformation of a subglacial till. Therefore, 
we infer that the base of Ruth Glacier, in the area in which the pulses occur, must be temperate. 
Additionally, theory predicts that the manner in which the glacier moved during the pulses, with 
a 20 km section of ice accelerating and decelerating in unison, and the lack of a wave front, is 
suggestive of a perennial temperate base (Fowler and others, 2001). The oscillating dynamic 
behavior is theorized to be due to a thin till with low permeability that is just above a critical 
value (Fowler and others, 2001). The basal roughness is theorized to be great enough to restrict 
ice flow and prevent the occurrence of major surges. These findings illustrate the influence of 
bedrock lithology and a glacier’s basal thermal regime on its dynamics, the combination of 
which will ultimately affect its mass balance. 
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545 


546 

547 Figure 1 . View of Ruth Glacier in the Alaska Range, Alaska. Red polygons indicate the extent 

548 of the granitic Tertiary age bedrock formations in the area, labeled Tpgr. KJf indicates 

549 Cretaceous/Jurassic age sedimentary bedrock, and Tty indicates Tertiary age sedimentary 

550 bedrock. The green line indicates a transect along the glacier centerline, beginning at the 

551 headwall of the north fork, measured in km. The background image was acquired by Landsat-5 

552 TM on Sept. 9, 1994. The black rectangle within the inset indicates the position of the 

553 background image within Alaska. 
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Figure 2. Velocity fields for Ruth Glacier, below the Great Gorge, from 1973 to 1983. Distances 
along the centerline transect are given in kilometers from the headwall of the north fork (See Fig. 
1). A) 1973-1974. B) 1974-1975. C) 1976-1977. D) 1977-1978. E) 1978-1980. F) 1980-1981. G) 
1981-1982. H) 1982-1983. 



Figure 3. Velocity fields for Ruth Glacier, below the Great Gorge, from 1983 to 1999. Distances 
along the centerline transect are given in kilometers from the headwall of the north fork (See Fig. 
1). A) 1983-1984. B) 1984-1985. C) 1985-1986. D) 1986-1987. E) 1987-1991. F) 1991-1992. G) 
1994-1995. H) 1995-1999. 
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570 Figure 4. Velocity fields for Ruth Glacier, below the Great Gorge, from 1999 to 2006. Distances 

571 along the centerline transect are given in kilometers from the headwall of the north fork (See Fig. 

572 1). A) 1999-2000. B) 2000-2001. C) 2001-2002. D) 2002-2003. E) 2003-2004. F) 2004-2005. G) 

573 2005-2006. 
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575 Figure 5. Velocity fields for Ruth Glacier, below the Great Gorge, from 2006 to 2012. Distances 

576 along the centerline transect are given in kilometers from the headwall of the north fork (See Fig. 

577 1). A) 2006-2007. B) 2007-2008. C) 2008-2009. D) 2009-2010. E) 2010-2011. F) 2011-2012. 
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Figure 6. Temporal evolution of ice surface velocity on Ruth Glacier at the 36 km and 43 km 
marks from 1973 to 2012. Five pulses are evident, with peaks in velocity in 1981, 1989, 1997, 
2003, and 2010, resulting in a characteristic pulse frequency of approximately 7 years. 
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Figure 7. Longitudinal profiles of ice velocity, acceleration, and strain rate along the centerline 
transects shown in Figure 5, for the years 2006 to 2012, which surround the 2009-2010 pulse. A) 
Ice velocity. B) Ice acceleration. C) Longitudinal strain rate. 
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602 Figure 8. Longitudinal profiles of ice velocity, acceleration, and strain rate along the centerline 

603 transects shown in Figure 4, for the years 1999 to 2006, which surround the 2002-2003 pulse. A) 

604 Ice velocity. B) Ice acceleration. C) Longitudinal strain rate. 
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Figure 9. Longitudinal profiles of ice velocity, acceleration, and strain rate along the centerline 
transects shown in Figure 3, for the years 1983 to 1999, which surround the 1987-91 and 1995- 
1999 pulses. A) Ice velocity. B) Ice acceleration. C) Longitudinal strain rate. 
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Figure 10. Longitudinal profiles of ice velocity, acceleration, and strain rate along the centerline 
transects shown in Figure 2, for the years 1973 to 1983, which surround the 1980-1981 pulse. A) 
Ice velocity. B) Ice acceleration. C) Longitudinal strain rate. 
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617 Table 1 . Landsat images used in this study. Upon receipt from the U.S. Geological Survey, Multispectral Scanner 

618 System images (MSS) have a spatial resolution of 60 meters, and the Thematic Mapper (TM) and Enhanced 

619 Thematic Mapper Plus (ETM+) images have spatial resolutions of 30 meters. 

620 


Image date 

Sensor 

Spatial resolution 
(meters) 

24 Sept. 1973 

MSS 1 

60 

27 July 1974 

MSS 1 

60 

23 Sept. 1975 

MSS 2 

60 

7 July 1976 

MSS 2 

60 

lJuly 1977 

MSS 2 

60 

2 Aug. 1978 

MSS 2 

60 

22 July 1980 

MSS 2 

60 

4 Aug. 1981 

MSS 2 

60 

3 July 1982 

MSS 3 

60 

18 Aug. 1983 

MSS 4 

60 

28 Aug. 1984 

MSS 5 

60 

16 Sept. 1985 

MSS 5 

60 

16 Sept. 1985 

TM 5 

30 

17 July 1986 

TM 5 

30 

21 Aug. 1987 

TM 5 

30 

29 June 1991 

TM 5 

30 

26 Aug. 1992 

TM 4 

30 

9 Sept. 1994 

TM 5 

30 

8 June 1995 

TM 5 

30 

19 June 1999 

TM 5 

30 

16 Aug. 2000 

ETM+7 

30 

2 July 2001 

ETM+7 

30 

21 July 2002 

ETM+7 

30 

9 Aug. 2003 

ETM+7 

30 

10 July 2004 

ETM+7 

30 

14 Aug. 2005 

ETM+7 

30 

10 Sept. 2006 

TM 5 

30 

28 Aug. 2007 

TM 5 

30 

5 July 2008 

ETM+7 

30 

8 July 2009 

ETM+7 

30 

20 Aug. 2010 

TM 5 

30 

6 July 2011 

TM 5 

30 

18 Oct. 2011 

ETM+7 

30 

20 Oct. 2012 

ETM+7 

30 
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629 Table 2. Georeferencing error and residual error for each image pair used to produce velocity fields. Georeferencing 

630 error represents the initial misalignment between the two satellite images as measured using COSI-Corr, and the 

63 1 residual error represents their misalignment after the georeferencing error is removed. Image geolocation RMSE is 

632 the root mean square error of the geometric misalignment of each image with respect to ground control points, as 

633 determined during processing of the images before distribution to the public domain. Values given are mean ± 1 

634 standard deviation. 


Image pair 

Sensor 

Georef. Error 

Georef. Error 

Residual Error 

Residual Error 

Image Geolocation 
RMSE (pixels, 
image year) 

types 

(pixels) 

(meters) 

(pixels) 

(meters) 

24 Sept. 1973, 27 July 1974 

MSS 

0.5010.21 

36.81 15.4 

0.26 + 0.14 

19.4 1 10.4 

0.386, 1973 

27 July 1974, 23 Sept. 1975 

MSS 

0.4410.15 

23.118.1 

0.20 + 0.12 

10.9 + 6.4 

0.396, 1974 

23 Sept. 1975, 7 July 1976 

MSS 

0.4110.61 

14.7121.7 

0.35 + 0.59 

12.5 + 20.9 

0.397, 1975 

7 July 1976, lJuly 1977 

MSS 

0.3810.10 

22.8 + 6.2 

0.14 + 0.07 

8.814.7 

0.428, 1976 

lJuly 1977, 2 Aug. 1978 

MSS 

0.33 10.13 

18.5 + 7.3 

0.16 + 0.07 

9.014.2 

0.391, 1977 

2 Aug. 1978, 22 July 1980 

MSS 

0.5010.24 

30.2 + 14.6 

0.26 + 0.12 

15.817.2 

0.357, 1978 

22 July 1980, 4 Aug. 1981 

MSS 

0.5010.22 

30.5 1 13.2 

0.26 + 0.15 

15.7 + 9.4 

0.682, 1980 

2 Aug. 1981, 3 July 1982 

MSS 

0.8710.87 

28.4 + 28.5 

0.36 + 0.34 

21.5 121.2 

0.747, 1981 

3 July 1982, 18 Aug. 1983 

MSS 

0.4210.19 

23.2 + 10.4 

0.32 + 0.17 

17.419.6 

0.430, 1982 

18 Aug. 1983, 28 Aug. 1984 

MSS 

0.2610.12 

15.8 + 7.3 

0.16 + 0.09 

10.115.9 

0.564, 1983 

28 Aug. 1984, 16 Sept. 1985 

MSS 

0.35 + 0.16 

21.1 + 10.0 

0.24 + 0.18 

14.6 + 11.1 

0.496, 1984, 
0.470, 1985 

Average for MSS imagery 


0.45 10.15 

24.116.7 

0.25 + 0.07 

14.1 14.2 

0.4710.12 

16 Sept. 1985, 17 July 1986 

TM 

0.45 10.25 

16.5 + 9.1 

0.38 + 0.24 

13.9 + 8.8 

0.228, 1985 

17 July 1986, 21 Aug. 1987 

TM 

0.4410.13 

12.2 + 3.7 

0.21 + 0.15 

5.714.3 

0.184, 1986 

21 Aug. 1987, 29 June 1991 

TM 

0.2010.14 

6.214.4 

0.16 + 0.09 

5.01 2.8 

0.155, 1987 

29 June 1991, 26 Aug. 1992 

TM 

0.2810.18 

7.815.0 

0.2110.13 

5.913.5 

0.191, 1991, 
0.132, 1992 

9 Sept. 1994, 8 June 1995 

TM 

0.2410.16 

8.816.1 

0.19 + 0.13 

7.115.0 

0.143, 1994 

8 June 1995, 19 June 1999 

TM 

0.2110.09 

6.5 1 2.8 

0.13 + 0.03 

4.01 1.1 

0.154, 1995 

19 June 1999, 16 Aug. 2000 

TM/ETM+ 

0.2210.16 

6.014.5 

0.15 + 0.13 

4.313.7 

0.155, 1999 

16 Aug. 2000, 2 July 2001 

TM 

0.2010.14 

7.014.9 

0.17 + 0.15 

5.815.1 

0.119, 2000 

2 July 2001, 21 July 2002 

TM 

0.3010.09 

9.01 2.7 

0.06 1 0.04 

1.911.2 

0.141, 2001 

21 July 2002, 9 Aug. 2003 

TM 

0.1410.10 

4.01 2.9 

0.12 + 0.11 

3.61 3.3 

0.187, 2002 

9 Aug. 2003, 10 July 2004 

TM 

0.1310.06 

3.91 1.8 

0.11 + 0.04 

3.31 1.4 

0.125, 2003 

10 July 2004, 14 Aug. 2005 

TM 

0.1310.07 

4.01 2.2 

0.10 + 0.07 

3.112.1 

0.151, 2004 

14 Aug. 2005, 10 Sept. 2006 

ETM+/TM 

0.2610.19 

7.215.3 

0.24 + 0.14 

6.61 3.9 

0.121, 2005 

10 Sept. 2006, 28 Aug. 2007 

TM 

0.1710.16 

5.3 1 5.0 

0.15 + 0.12 

4.61 3.6 

0.165, 2006 

28 Aug. 2007, 5 July 2008 

TM/ETM+ 

0.1810.08 

6.512.8 

0.16 + 0.09 

5.713.4 

0.127, 2007 

5 July 2008, 8 July 2009 

TM 

0.1610.07 

4.812.1 

0.06 + 0.03 

2.010.9 

0.148, 2008 

8 July 2009, 20 Aug. 2010 

ETM+/TM 

0.1810.06 

5.61 1.9 

0.09 1 0.04 

2.71 1.4 

0.139, 2009 

20 Aug. 2010, 6 July 2011 

TM 

0.1810.13 

5.614.1 

0.1110.11 

3.513.4 

0.092, 2010, 
0.225, 2011 

18 Oct. 2011, 20 Oct. 2012 

ETM+ 

0.2110.14 

6.3 14.3 

0.18 + 0.14 

5.514.4 

0.243, 2011, 
0.290, 2012 

Average for TM/ETM+ imagery 


0.2210.09 

7.01 3.0 

0.15 + 0.07 

4.912.6 

0.16 ± 0.04 


32 


635 


33 


